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Abstract 

We present a lattice computation of the first moment of the kaon's leading- 
twist distribution amplitude. The results were computed using ensembles 
with 2+1 dynamical flavours with the domain wall fermion action and 
Iwasaki gauge action from the RBC and UKQCD joint dataset. The first 
moment is non-zero because of S'f/(3)-breaking effects, and we find that we 
are able to measure these effects very clearly. We observe the expected chi- 
ral behaviour and finally obtain (f )(2GeV) = 3/5a^(2GeV) = 0.032(3), 
which agrees very well with results obtained using sum-rules, but with a 
significantly smaller error. 
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1 Introduction 

Hadronic light-cone distribution amplitudes are fundamental non-perturbative ingredients 
in the QCD analysis of hard exclusive processes. Phenomenological applications which 
require knowledge of the distribution amplitudes include electromagnetic form-factors 
at large momentum transfer and related processes [1-7]. More recently, following the 
development of the factorization framework, the distribution amplitudes are also an im- 
portant component in the phenomenology of exclusive charmless two-body 5-decays (i.e. 
.B-decays into two light mesons) [8-14]. These are a particularly important set of pro- 
cesses for CKM-analyses and for studies of CP-violation. Here, we present the first results 
from our new lattice project in which we are computing the moments of the light-cone 
distribution amplitudes for light pseudoscalar and vector mesons. 

The subject of this letter is the first moment of the leading-twist distribution amplitude 
of the kaon, 4>k{u, fi), which parametrizes the overlap of a kaon with longitudinal momen- 
tum p with the lowest Fock state consisting of a quark and an anti-quark carrying the 
momentum fractions up and up — (1 — u)p, respectively (u + u — 1). It is defined by the 
non-local (light-cone) matrix element: 



represents the path-ordered exponential from —z to z, so that the bi-local current on the 
left-hand side of eq.(l) is gauge invariant. The distribution amplitude is normalized by 
Jg 1 du 4>k(u, //) = 1, and can be expanded in terms of Gegenbauer polynomials Cn 2 (2u— 1), 



The lowest-order anomalous dimensions of the moments a^(fi) grow with n and thus 
higher moments may be suppressed for large values of the renormalization scale /i. In this 
paper we present our results for the lowest Gegenbauer moment af - , which is proportional 
to the average difference of the longitudinal quark and anti-quark momenta of the lowest 
Fock state: 



(0\q(z)^ 5 V(z,-z)s(-z)\K(p))\ z2=Q = f K (ip p ) [ due^- u > z Mu,^), (1) 









(4) 



af = 5/3(£) is obtained from the matrix element of a local operator, 
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eq.(5) is the leading term in the Taylor expansion of expression (1) around z — 0. Our 
conventions for the covariant derivatives are = _D M — D^, = d ^ + igA^ and 
0^ = 3^- igA^. 

The first moment of the kaon's distribution amplitude has in the past been determined 
mainly from QCD sum rules, and recent results include: 

af(lGeV) = 0.05(2) [15], 0.10(12) [16], 0.050(25) [17] and 0.06(3) [18] . (6) 

In this work we present the results of a lattice study of this quantity using Nf = 2 + 1 
dynamical flavours of domain wall fermions [19,20], which have good chiral properties, 
and the Iwasaki gauge action [21,22]. We used gauge field ensembles from the RBC and 
UKQCD dataset with three values of the light-quark mass and the calculations are carried 
out with equal valence and sea quark masses (i.e. with full unitarity). Further details of 
the simulation can be found in sec. 4 below. We have observed a clear signal for partonic 
S'6 r (3)-breaking effects in the leading-twist kaon distribution amplitude and find that the 
first moment satisfies the chiral behaviour expected from chiral perturbation theory. For 
our best results we quote 

(£) M %u = 2GeV) = 0.032 ±0.003 (af (2 GeV) = 0.053(5)) (7) 
(£) m (li = lGeV) = 0.040 ±0.004 (af (1 GeV) = 0.066(6)) , (8) 

in agreement with most previous results. The errors are already small, and, as discussed 
below, will decrease still further in the near future. 

While we were completing this paper, Braun et al. released the following result from a 
lattice simulation using improved Wilson fermions [23]: 

af (2 GeV) = 0.0453 ± 0.0009 ± 0.0029 , (9) 

in reasonable agreement with our result. 

The plan of the remainder of this paper is as follows. In the following section we introduce 
the basic definitions, in particular the ratio of Euclidean lattice correlation functions from 
which we determine ( £ ) and which we estimate using a Monte Carlo simulation. Section 3 
contains the discussion of the perturbative calculation of the renormalization constants. 
We present our results for the bare matrix elements in section 4 and combine them with 
the renormalization constants and discuss the systematic uncertainties in section 5, where 
we present our final result. We end with a brief summary and conclusions (sec. 6). 



2 ^£^bare f rom Lattice Correlation Functions 

We start by briefly describing the overall strategy of our calculation of ( £ ). As we explain 
in this section, we exploit the fact that ( £ ) can be obtained directly from a ratio of two 
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Euclidean correlation functions, which we evaluate using a Monte Carlo simulation. The 
statistical fluctuations are reduced in the ratio and, for each choice of quark masses used 
in the simulation, we are able to obtain ( £ ) with good precision. 

In constructing the lattice operators, we use the following symmetric left- and right-acting 
covariant derivatives: 

D^(x) = ^- { U(x, x + p,)i/>(x + p) -U(x,x - fi)ip{x - p)} , (10) 

and ^ 

ip(x)D tM = — + p)U(x + ji, x) — ip(x — fi)U{x — £i,x)} . (11) 

where the C/'s are the gauge links and fi is a vector of length a in the direction fi (a 
denotes the lattice spacing). 

To illustrate the method, consider the local lattice operators 1 : 

<— > 

O pix (x) = g(a;)7 P 75 s(x) , A p (x) = q(x) 7 P 7 5 s(x) and P(x) = q(x) 75 s(x) , (12) 
and define the two-point correlation functions 

C P »(t,p) = J2e^(0\O p ,(t,x )Pt(0)|0> , (13) 

X 

and 

C Av p(t,p) = J2e^(0\A v (t,x )P t (0)|0) . (14) 

X 

Here q and s represent the light and strange quark fields, respectively. At large Euclidean 
times t and T — t, where T is the length of the lattice in the time direction, the correlation 
functions (13) and (14) approach 

hi 



CA p{t p) ^ z pf ^ e -^nL m _ TmEK) x fe) 

where we have used (4), (5) and defined Z P = (i\"(p)|Pf |0) and the kaon's bare decay 
constant (0\A I/ \K(p)) = ip^f^ . The superscript frare denotes the fact that the operators 
are the bare ones in the lattice theory with ultraviolet cut-off a -1 . Taking the ratio of the 
two correlation functions 

\ bare / 1 rj\ 

^A„P{t,P) Pu 

allows us readily to extract the bare value of (£). The braces in the subscripts {pfu,} 
indicate that the indices are symmetrized. 

1 In the simulation we actually use a smeared pseudoscalar density P in order to improve the overlap 
with the kaon state. The present discussion holds for both local and smeared pseudoscalar densities. 
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Figure 1: One-loop vertex diagrams evaluated in the perturbative renormalization of 0{ pp j 




Figure 2: One-loop diagrams contributing to the quarks' wave function renormalization. 



In the continuum 0{ pp } transforms as a second rank Lorentz tensor, whereas on the 
lattice we need to consider transformation properties under the hypercubic group and 
discrete symmetries, in particular parity and charge conjugation. We choose to evaluate 

the matrix elements of 0{ pp } for /i 7^ p. These operators transform as a six-dimensional 
representation of the hypercubic group and lattice symmetries exclude mixing with other 
operators [24] including operators containing a total derivative (since we are considering 
matrix elements with a non-zero momentum transfer, operators with a total derivative 
have to be considered). As can be seen from eq. (17), in order for the matrix elements 
to be non-zero, we require both p p 7^ and p p 7^ 0. We satisfy this condition by taking 
/i — v — 4, p — 1, 2 or 3 and \p p \ = 2n/L . 

Having obtained (£) bare > we need to determine the renormalization factor relating the 

lattice operators and A v to the corresponding continuum operators in some standard 
renormalization scheme; here we do this at one-loop order in perturbation theory as 
explained in the following section. 

3 Perturbative Renormalization of the Lattice Oper- 
ators 

The perturbative matching from the lattice to MS schemes is performed by comparing 
one-loop calculations of the amputated two-point Green function with an insertion of 
the operator 0{ pp j in both schemes (which requires the evaluation of the diagrams in 
fig. 1), together with appropriate wave function renormalization factors (fig. 2). Defining 
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0|^}(/i) = Z 0{pM} 0^"|(a), the renormalization factor is given by 



Zo { 



«*> (1 - < 



1 + ("I ln(/i2fl2) + " Sl + ^ " V 



(18) 



In this expression, (1 — Wq)Z w is a characteristic normalization factor for the physical 
quark fields in the domain wall formalism. It is a common factor in the numerator 
and denominator of the ratio R{ Pil y tV as are the contributions from the wave function 
renormalization. Z w represents an additive renormalization of the large Dirac mass or 
domain wall height M — 1 — w which can be rewritten in multiplicative form at one-loop 
as 

7 - i i 9 2 C F 2w 

-L ~r ^ w 1 -...2 w ' V / 



^0 

The one-loop correction z w becomes very large for certain choices of M [25,26], including 
that used in our numerical simulations, so that some form of mean-field improvement is 
necessary, as discussed below 2 . 

The terms Ep and Si come from quark wave function renormalization. The terms V MS 
and V come from the one-loop corrections to the amputated two-point function. They 
are given by "vertex" and "sail" diagrams, plus an operator tadpole diagram in the lattice 
case. Using naive dimensional regularization (NDR) in Feynman gauge with a gluon mass 
infrared (IR) regulator, 

e ms = 1 ? l^=-|. (20) 

The lattice contribution Si has been evaluated for domain wall fermions with the Iwasaki 
gluon action in Feynman gauge and a gluon mass IR regulator in [26]. We have calculated 
the lattice vertex term V for the same action, gauge and IR regulator to complete the 
evaluation of Zq,, • The perturbative calculation is explained in [25-27] and the form of 
the Iwasaki gluon propagator can be found in [28]. Values for V are given as a function 
of M in Table 1. Chiral symmetry of the domain wall action implies that these results 
also apply for the operator which is 0{ PM } without the 75. We have confirmed that our 
results reproduce those found by Capitani [27] if we replace the gluon propagator for 
the Iwasaki gauge action by the propagator corresponding to the standard plaquette 
action. This provides a powerful check of our calculation. We note that the perturbative 
renormalization factor for the same operator using overlap fermions and the Liischer- 
Weisz gauge action can be found in [29]. 

Our numerical simulations use M = 1.8. For this value of M, with the Iwasaki gluon 
action, the one-loop coefficient z w in the physical quark normalization can be extracted 
from S w in Table III of [26]. This gives z w 112, making it clear that mean-field 
improvement is necessary. We follow the prescription described in [26]. 

2 The factor 1/(1 — Wq)Z w cancels however, in the evaluation of the ratio Zo {p ^ } /Za, where Za is the 
renormalization constant for the axial current. 
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4.0343 
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-0.8050 
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-0.9135 
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1.8 


111.69 


37.901 


4.3249 


-3.9462 


-1.2676 


1.9 


241.55 


84.270 


4.5209 


-3.8447 


-1.5651 



Table 1: Constants needed for the perturbative renormalization of the operator 0{ p ^ using 
domain wall fermions and the Iwasaki gauge action (ci = —0.331). M is the domain wall 
height and #£1 + SV = E^ ~ ^i + ^ MS — V ■> while other quantities are defined in the 
text. Si and V are dependent on the gauge and the infrared regulator: Feynman gauge 
and a gluon mass are used here. zi MF ' ) and Ei are extracted from the results in [26], while 
V has been calculated as part of this work. 



The first step is to define a mean- field value for the domain wall height, 

m mf = jvf — 4(1 — P 1/4 ) (21) 
where P = 0.58813(4) is the average plaquette value in our simulations. This leads to 

m mf = ! 3029. (22) 

The physical quark normalization factor becomes [l — (wJf F ) 2 ] Z^ F , with 

= l + f^hr, C = l _ ^mf )2 + 3 2-^mf) = 5.2509, (23) 

where Tmf = 0.0525664 [26] is a mean-field tadpole factor and E w is evaluated at M MF , 
leading to z™ F = 5.2509. Likewise, E x = 3.9731 and V = -4.1907 in equation (18) are 
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evaluated at M MF and the mean-field improved renormalization factor for our simulations 
becomes: 



o 



0.9082 



(24) 



We make two choices for the mean-field improved MS coupling. The first uses the mea- 
sured plaquette value, P, according to [26] 

IP 22 

+ d g + c p + -— In (//a) , (25) 



W 9 2 9 P l&r* 

where d g = 0.1053 and c p = 0.1401 for the Iwasaki gauge action and (3 = Q/g 2 = 2.13 in 
our simulations. The second choice is the usual continuum MS coupling. At fia = 1, we 
find aj^g(plaq) = 0.1752 and a^jg(ctm) = 0.3385. This disparity in the values of the two 
couplings is further motivation for the programme of non-perturbative renormalization 
which we are currently undertaking. With these two choices of coupling, our value for the 
renormalization factor becomes: 

J 0.9811 plaquette coupling 
° {PM} ~ [0.8719 continuum MS. ^ ' 

We also evaluate the mean-field improved expression for the axial vector current [26], 
interpolating to our mean- field M MF , and obtain 

0.7947 plaquette coupling 
0.6514 continuum MS . 

The ratio of the two renormalization factors is 

-Zo {pm} J 1.2346 plaquette coupling 
1.3384 continuum MS . 

For the purposes of this letter we include the spread of results in eq.(28) as the estimate 
of our current systematic uncertainty of the renormalization factor 3 . This uncertainty 
will be significantly reduced as we complete our programme of non-perturbative renor- 
malization. (£) bare should be multiplied by the factor on the right-hand side of eq.(28) 
to obtain the result in the MS scheme. For this factor we take 

= 1-28 ±0.05. (29) 



3 We mention in passing that using the bare coupling, = 1-26 • 
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4 Numerical Simulation and Results 



The numerical results presented here are based on Monte Carlo estimates of correlation 
functions evaluated on representative sets of UKQCD/RBC gauge field configurations 
that were generated with Nf — 2 + 1 flavours of dynamical domain wall fermions [19, 20] 
with Iwasaki Gauge action [21,22] using the QCDOC computer [30-33]. The hadronic 
spectrum and other properties of these configurations have been studied in detail and the 
results will be presented in ref. [34]. We wish to thank our collaborators in the RBC and 
UKQCD collaborations for access to several preliminary results which we require for this 
study. These are specifically m res in the chiral limit, the bare kaon pseudoscalar masses, 
and the inverse lattice spacing. 

The preliminary nature of these intermediate values does not introduce a significant un- 
certainty in our results compared to our statistical errors and the systematic errors due 
to the chiral extrapolation and renormalization. Indeed we shall only quote the lattice 
spacing to two significant figures and demonstrate the insensitivity of our results to this. 
However, we caution the reader that definitive results for these intermediate quantities 
will be found in the forthcoming publication. 

The lattice volume is (L/a) 3 x T/a = 16 3 x 32 and the length of the fifth dimension 
is L s = 16. The choice of bare parameters is f3 — 2.13 for the bare gauge coupling, 
am s = 0.04 for the strange quark mass (which has been tuned to correspond to the 
physical value) and am q = 0.03, 0.02, 0.01 for the bare light-quark masses. With this 
choice of simulation parameters the lattice spacing is cr 1 = 1.6 GeV [34]. Due to the 
remnant chiral symmetry breaking the quark mass has to be corrected additively by the 
residual mass in the chiral limit, am res = 0.003 [34]. 

Statistical errors for observables have been estimated both with the jack-knife and with 
a direct estimation of the integrated autocorrelation time as suggested in [35]. 

4.1 Bare correlation functions 

For each value of the light-quark mass we computed the correlation functions on 300 
gauge configurations separated by 10 trajectories in the Monte Carlo history. On each 
configuration we averaged the results obtained from 4 sources for the lightest quark mass 
(m q a = 0.01) and 2 sources for the remaining two masses (m q a = 0.02 and 0.03). The 
sources were chosen to be at the origin and (8,8,8,16) for all three masses, and in addition 
at (4,4,4,8) and (12,12,12,24) for m q a = 0.01. In order to improve the overlap with the 
ground state, at the source where we insert the density P\ we employed gauge invariant 
Jacobi smearing [36] (radius 4 and 40 iterations) with APE-smeared links in the covariant 
Laplacian operator (4 steps and smearing factor 2) [37,38]. 

The preliminary kaon masses corresponding to the simulated bare light-quark masses are 
am om = 0.4164(10), am° K 02 = 0.3854(10), and am° K 01 = 0.3549(14) [34]. 
In order to extract (£) from the ratio R{ pfl y. u (t,p) defined in (17) we need \p\ ^ 0. 
Since hadronic observables with larger lattice momenta have larger lattice artefacts and 
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Figure 3: Jack-knife results for (£ ) bare as a function of the time. The ranges over which 
we fit and the corresponding results are indicated by the black lines. 



am ud 



0.03 



0.02 



0.01 



(u-u) 0.0057(4) 0.0119(10) 0.0181(18) 



X- limit 



0.0262(23) 



Table 2: Summary of results for the bare values of the 1st moment of the kaon's distri- 
bution amplitude. The result we obtain after the linear chiral extrapolation is quoted in 
the right-most column. 



statistical errors, we restrict the choice of indices to p = v = 4 4 and \i = 1,2,3 with 
\p | = 2tt I L = ±2tt/L with the remaining two components of p equal to 0). ( £ ) bare is 
then obtained from the correlation function at large times: 

R m -At, Pk = ±2tt/L) = ± i — ( £ ) barc , \p | = — , k = 1, 2, 3 . (30) 

The plot in figure 3 shows our results for ( £ ) as a function of t obtained from the ratio 
R{4k}-,4{t,pk = ±2ii/L) for the three values of the mass of the light quark. The results 
have been averaged over the three values for k and the 6 equivalent lattice momenta with 
\p\ = 2tt/L and combining the results at t with those at T — t — 1. There are clear 
plateaus, demonstrating that the S"t/(3)-breaking effects are measurable and ( £ ) can be 
determined. The results from the fits for (£ ) bare are summarized in table 2. 

4 The index 4 corresponds to the time-direction. 
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4.2 Chiral extrapolation 



For the pion the first moment a\ vanishes since isospin symmetry induces invariance 
under u <-> (1 — u). For a non-degenerate quark- ant iquark pair (such as the kaon) flavour 
symmetry breaking implies that the first moment of the distribution amplitude is non- 
zero. The leading SU (3)- violating effects for the kaon's distribution amplitude have been 
studied at next-to-leading order in chiral perturbation theory (%PT) in [39,40]. No chiral 
logarithms appear and the prediction for the mass-dependence is 

(0 = ~jr( ms - m ud)bi,2, (31) 

where / and B are conventional %PT parameters and &i )2 is a Wilson coefficient as 
introduced in [39]. 

We plot our results for (£ ) bare as a function of the light-quark mass in fig. 4. We take 
into account the remnant chiral symmetry breaking by defining the chiral limit at the 
point am q + am rcs = 0. The linear behaviour predicted in eq. (31) is well satisfied (with 
a tiny x 2 /d.o.f. of about 10" 5 ). Moreover the line passes through (£ ) barc = at a value 
of the light-quark mass (denoted by the open square in fig. 4) which is consistent with 
the mass of the strange quark, as expected for the 577(3) symmetric case {am U( i = am s = 
0.04). More specifically, the intercept of the linear fit with the (£) bare = ax i s occurs at 
am ud = 0.0391+^ 3 7 . 

From the linear fit 5 we obtain (£) bare = 0.0262(23) in the chiral limit and in the next 
section we combine this result with the renormalization constant in eq.(29) to arrive at 
our final result. 



5 Systematic Uncertainties and our Final Result 

To obtain the final result for the (£) in the MS scheme at // ~ 1.6 GeV we multiply the 
result obtained from the bare operators with cut-off a~ l = 1.6 GeV, (£ ) bare = 0.0262(23), 

by the ratio of renormalization factors in eq.(29), ^ = 1.28(5): 

(£) m (n = 1.6 GeV) = 0.034 ±0.003. (32) 

In order to be able to compare our result with previous calculations we evolve it to 
renormalization scales of 1 GeV and 2 GeV using the three-loop anomalous dimension [41], 
obtaining: 

(£) H %it = 2GeV) = 0.032 ±0.003 (33) 

(Z) m (n = lGeV) = 0.040 ±0.004. (34) 

5 We have also performed quadratic fits to the chiral behaviour but the results do not change in any 
significant way. 
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The error in the renormalization factor due to the uncertainty in the lattice spacing is 
negligible. For example if we conservatively allow the lattice spacing to vary between 
1.58 GeV and 1.62 GeV, the contribution to the relative error on (£) is less than 0.2%. 

Among the uncertainties which, at this stage at least, we are not in a position to check 
numerically are the continuum extrapolation, finite-volume effects and the fact that the 
strange quark mass (m s a = 0.04) is only approximately tuned to its physical value. The 
lattice artefacts are formally of 0{a 2 K^ CD ) ~ 2.5%. We would expect the finite volume 
effects to be small and are currently checking this with a simulation on a 24 3 x 64 lattice. 
The strange quark mass appears to be well tuned [34] so again we expect the contribution 
to the error from this uncertainty to be very small. Thus we expect the errors from these 
three sources to be sufficiently small not to change the errors quoted in eqs.(33) and (34) 
which we take to be our best estimates. We are also carrying out a systematic programme 
of non-perturbative renormalization which will enable us to reduce the uncertainty in the 
renormalization constants. 



6 Summary and Conclusions 

In this letter we have presented the first results from our major lattice study of the 
leading-twist light-cone distribution amplitudes of the light mesons. We demonstrate 
that the SU (3)-breaking effects which lead to a non-zero value for the first moment of the 
kaon's distribution amplitude are sufficiently large to be calculable in lattice computations 
and satisfy the expected chiral behaviour. Our results for the first moment are presented 
in eqs.(33) and (34). 
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We are in the process of implementing a number of improvements, including non-pert- 
urbative renormalization and a simulation on a larger lattice (24 3 x 64) which will help 
reduce these two sources of systematic error (or in the latter case give us confidence in our 
expectation that the finite- volume errors are indeed small). We will also produce results 
for the second moment of the pion's and kaon's distribution amplitudes and for those of 
the p and K* vector mesons. In order to reduce the lattice artefacts we will investigate 
the use of partially twisted boundary conditions [42-44] which will allow us to calculate 
the observables at smaller values of lattice momenta. 

While we were completing this paper, ref. [23] appeared with the results of an Nf = 
2 study of the distribution amplitude, using 0(a) improved Wilson fermions and the 
plaquette Wilson gauge action. For the first moment of the kaon's distribution amplitude 
the calculation was performed at f3 — 5.29 corresponding to a lattice spacing of about 
2.6 GeV. The result of ref. [23] is given in eq. (9) above. 
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